Detection of an X-ray Pulsation for the Gamma-ray Pulsar 

Centered in CTA 1 
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ABSTRACT 

We report the detection of X-ray pulsations with a period of ^315.87 ms from the 2009 
XMM-Newton observation for the radio-quiet 7-ray pulsar, LAT PSR J0007+7303, centered in 
the supernova remnant CTA 1. The detected pulsed period is consistent with the 7-ray periodicity 
at the same epoch found with the Fermi Gamma-ray Space Telescope. The broader sinusoidal 
structure in the folded light curve of the X-ray emission is dissimilar to that of the 7-ray emission, 
and the phase of the peak is about 0.5 shifting from the peak in the 7-ray bands, indicating that 
the main component of the X-rays originates from different sites of the pulsar. Wc conclude 
that the main component of the X-ray pulsation is contributed by the thermal emission from 
the neutron star. Although with a significantly different characteristic age, PSR J0007+7303 is 
similar to Geminga in emission properties of X-rays and 7-rays; this makes PSR J0007+7303 the 
second radio-quiet 7-ray pulsar with detected X-ray pulsations after Geminga. 



Subject headings: radiation mechanisms: 
individual (PSR J0007+7303, Geminga) - 

Introduction 



non-thermal — radiation mechanisms: thermal — pulsars: 
- X-rays: general — gamma rays: general 



In the first year of its operation, Fermi ob- 
servatory detected 23 new radio-quiet 7-ray pul- 



sars (|Abdo et al.ll2009HSaz Parkinson et alll2010l) 

The X-ray counterparts for some of these 7-ray 
pulsars have been sugg estively proposed, e.g ., 



and 2XMM J202131. 0+402645 (|Trepl et al 



CXOU J180950.2-23 3223 dBraie fc Romani 
RX J1836.2+5925 (Irlalpern et all 120021 



2002), 



2007), 



2010) 



these sources show similar spectral features as 
a neutron star. Although these p ulsars might 



be cl aimed as a next Geminga (e.g., lAbdo et al 
20101) . their pulsed emission in X-ray bands have 
not been confirmed so far. 
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The X-rays from these radio-quiet 7-ray pul- 
sars are expected to be generally characterized 
by thermal emission from the neutron star sur- 
face and the hot polar cap as well as non-thermal 
emission from the magnetosphere and th e pulsar 
wind nebula (|Caraveo et al.ll2003l 120041) . These 
different emission components can possibly have 
thei r contributions in d ifferent rotational phases 
(cf. ICaraveo et al. 2004 ). Therefore, it is not pos- 
sible to obtain a complete picture of this unique 
class of pulsars without the detection of X-ray pul- 
sations. Nevertheless, until now, Geminga is still 
the only one in this class with the X-ray pulsed 
signals identified. 

In the search for X-ray pulsations from 7- 
ray-only-emitting pulsars, the point source cen- 
tered in the composite supernova remnant (SNR) 
CTA 1 (G119.5+10.2) is one of the most in- 
triguing candidates. CTA 1 was firstly dis- 
covere d by Owens Valley Rad io Observatory in 
1960s (lHarris fc Roberts! 119601 ); previous studies 
(Se ward et al J ll995l) revealed that the X-rays from 
the center of CTA 1 consist of diffuse emission plus 
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Fig. 1. — L eft panel: Sp i n periods of PSR J0007+7303 from Fermi and XMM. The squares represent the spin periods 
detected by I Abdo et all d2008L lioogl). T he plus sign marks the glitch of PSR J0007+7303 with df/f = 5 x 10~ 7 
(jAbdo fc Fermi-LAT Collaboration! l201Cl ). The circles represent the spin periodicities detected by the Fermi data 
before the glitch, and the solid circles represent those after the glitch. The triangle represents the period detected by 
the XMM data of 2009 with an u ncertainty shown b y the error bar. The dashed line represents the period derivative 
-3.6133 x 10 12 Hz/s derived bv lAbdo et "all (|200Sh . The periods after glitch have a little shift from the dashed 
line which are clearly shown in the insert panel. Right panel: Rayleigh test of XMM data around the frequency 
[3.114,3.218] Hz. The width of each independent trial is about 1.037 x 10 -5 Hz. The most significant signal appears 
at the frequency of 3.165844 Hz. 



a faint point source RX J0007.0+7302. The pre- 
cise position of RX J0 007. 0+7302 was pro vided 



by the Chandra image (jHalpern et alJ 120041 ). No 



radio counterpart has bee n found at the pos i 



tion of RX J0007.0+7302 ( Pineault et alJ ll993) 



On the other hand, the high 7-ray-to- X-ray and 



X-ray-to-optical flux ratios ([Brazier et al.l 11998 ; 



Halpern et al.ll2004 ) and the morphology of a bent 
jet from a torus-like pulsar wind nebula (PWN) 



revealed in X-ray image (jHalpern et al.l 120041 ) 



all characterize RX J0007. 0+7302 as a neutron 
star. Further studies of the spectral properties 
for RX J0007.0+7302 were also consistent with 
a pulsar without radio emission, e.g., Gcminga 
( Slane et al. 20041 ). In addition, based on the po- 
sitional coincidence and the emission properties, 
RX J0007.0+7302 was also identified as the X- 
ray counterpart of 3EG J0010+7309, which is a 
7-ray source found by the Energetic Gamma-Ray 
Experiment Telescope (EGRE T) on the Corny 



ton Gamma-Raii O bservatory (ISlane et al.l 11997 



Brazier et aTlll998l) 



Although all these observations suggested that 
RX J0007.0+7302 is likely to be a radio-quiet 
7-ray pulsar associated with the supernova rem- 
n 19971 ). the final confirmation 
from the detection of 7 pulsations with 
The detected 7— ray 



nant (ISlane et al 
came 



Fermi (|Abdo et all 120081) 



pulsation has a spin period of 315.86 ms. This 
pulsar has a characteristic spin-down age of 
r ~ 14000 yr, which is consiste nt with the age 



of C TA 1 (10000 - 20000 yr, IZhang & Ghens 



1998 and ISlane et al.l 120041 ). This further sup- 



ports the scenario that this pulsar is indeed as- 
sociated with CTA 1. In addition, a precise po- 
sition of LAT PSR J0007+ 7303 can be inferred 
from the Chand ra image (jHalpern et al. 2004: 
Abdo et al. 20091 ). However, the search of the 



corresponding pulsed signal in the archival X- 
ray data obtai ned before 2009 was unsuccessful. 
Caraveol (|2009h indicated that the X-ray pulsation 
can be derived from the contemporaneous Fermi 
ephemeris of PSR J0007+7303 although no any 
further information of the X-ray pulsation was 
provided. This leaves the X-ray temporal proper- 
ties of RX J0007.0+7303 still unconstrained. 

In this Letter, we report the detection of X-ray 
pulsations from PSR J0007+7303 using the 2009 
XMM-Newton (hereafter XMM) observation (PI: 
Caraveol [20081 . ObsID: 06049401). We show de- 
tailed information of temporal and spectral analy- 
ses for this source. We investigate the source spec- 
trum and the folded light curves in different energy 
bands to determine the possible origins of the X- 
ray emission from PSR J0007+7303. 
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2. Observations and data analysis 

The central region of CTA 1 was observed by 
XMM twice with an epoch separation of ~ 7 
years. The first observation was performed on 
2002 February 21 for ~ 40 ks and the second one 
was carried out on 2009 March 7 for ~ 120 ks. 
These two observations were operated in the full- 
frame mode for the two MOS detectors and in 
the small-window mode for the pn detector. The 
medium and thin filters were used for the two MOS 
detectors and the pn detector, respectively. The 
temporal resolution of ^6 ms for the small- window 
mode of pn detector is suitable for the periodic- 
ity search o f X-ray pulses from PSR J0007+7303. 
Slane et al.l ([20041 ) have already analyzed the old 
XMM data observed on 2002 February 21 but no 
significant pulsation was detected using a total of 
1055 counts in the energy band 0.3-10 keV and an 
upper limit of 61% for the pulsed fraction was re- 
ported with a sinusoidal profile. In our study, we 
focus on the 2009 XMM data. All the data were 
processed with the SAS version of 10.0.0. 

2.1. Timing analysis 

We first performed the source detection by us- 
ing the task "edetect_chain" and determined the 
position of PSR J0007+7303 in the X-ray band at 
(J2000) R.A.=00 ft, 07 m 01 s .2, Dcc.=+73°03'07".6 
with an uncertainty of ~ 0".6. The point source 
was then extracted at this position within a circu- 
lar region of 15" in radius, which is consistent with 
the 70% of the encircled energy function. After 
screening the background flare, the effective expo- 
sure is ~ 67 ks. Following the standard process 
of the data reduction, 2989 counts were yielded 
in the energy band 0.2-12 keV. Furthermore, the 
photon arrival times were corrected to the solar 
system barycenter with the "barycen" task (JPL 
DE200 earth ephemeris) of the SAS package. 

lAbdo et al.l (|2009f ) reported a periodic signal of 
3.1658891845(5) Hz at the reference epoch MJD 
54754 with a frequency derivative of —3.6133(3) x 
10~ 12 Hz/s for PSR J0007+7303. Accord- 
ing to this timing ephemeris, a periodicity of 
3.165844344(4) Hz is expected at the start of 
the GTI of the 2009 XMM data (epoch MJD 
54897.6335129). We thus search for periodic sig- 
nals with the known frequency derivative around 
the predicted frequency in the X-ray data. By 




Fig. 2.— Folded light curves of PSR J0007+7303. 
The upper panel shows the pulse profile in the 0.2-12 
keV from the 2009 XMM observation and the lower 
one shows the profile in the 100 MeV to 300 GeV from 
the Fermi observations. The counts were normalized 
to the average photons of each bin. All the curves were 
folded simultaneously with the periodic frequency and 
the frequency derivative of the same epoch at the start 
of GTI of the 2009 XMM data as we described in the 
main text. 



performing an ff-test ( de Jager et al. 1989h . 



we 

found a test statistic of H ~ 27.0, which corre- 
sponds to a random probability of 2.0 x 10 -5 for a 
single trial ( de Jager fc Biisching 2010l ). The un- 
certainty of our detec tion was dete rmined by the 
equation provided by iLeahvl (|l987h and the most 
significant trial for the periodicity search was de- 
termined at 3.165844(5) Hz as shown in the left 
panel of Fig. [TJ The TJ-statistic also obtains the 
most significant periodic detection with a sinu- 
soidal structure consisting of only one harmonic. 
We then verified the Z\ value of 10000 indepen- 
dent nearby trials (~ 3.114 Hz - 3.218 Hz) by 
a Rayleigh test. Our detection of the periodic 
signal is highly significant as shown in the right 
panel of Fig. Q] 

In order to compare the pulse profiles in X-rays 
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and in 7-rays, we obtain ~ 3-month Fermi data 
(2009 January 21 - 2009 April 13) within the en- 
ergy range of 100 MeV to 300 GeV on the epochs 
close to the 2009 XMM observation. We selected 
a 1-degrce radius region centered on the source 
at R.A.=1.7565°, Dec.=73.0523° (J2000) accord- 
ing to the position of the X-ray counterpart de- 
tected by Chandra. The pulse profiles folded si- 
multaneously with the frequency of 3.165844344 
Hz and the frequency derivative of —3.6133(3) x 
10" 12 Hz/s at epoch MJD 54897.6335129 for the 
XMM and the Fermi observations are shown in 
Fig. [2] The pulse profile in the X-ray band shows 
a large pulsed phase (0.5625 - 1.3125) with 0.75 
duty cycle while in the 7-ray band the pulse profile 
shows a small pulsed phase (0.125 - 0.5625) with 
~ 0.45 duty cycle. In addition, the peaks of these 
two light curves shift about 0.5 phase, indicating 
that the X-rays and 7-rays originate from different 
sites of the pulsar. 

2.2. Spectral analysis 

The EPIC spectra of PSR J0007+7303 were 
extracted using the aforementioned 15"-radius 
source region. The background was defined as an 
annulus with the inner and outer radii of 75" and 
120" centered at the source region. The photon 
redistribution matrices and the ancillary region 
files were generated by the task "rmfgen" and 
"arfgen" of the SAS package. In order to compare 
with t h e previous analysis re ported bv lSlane et al 



Table 1: Best- fit spectral parameters for PSR 
J0007+7303. 



( 2004 ); Halpern et al.l ( 20041 ). our spectra were re 
stricted in the energy range 0.5 — 10.0 kcV and 
each spectral bins were rcbinned with a mini- 
mum of 30 counts. To reduce the uncertain- 
ties in spectral fitting, we set the column den- 
sity to be ~ 2.8 x 10 21 cm~ 2 according to the 
previous ASCA/R OSAT measurements of CTA 1 



(plane et al.l 119971 ). This v alue is also consisten t 



with the optical extinction (IHalpern et al.ll2004l) . 

We have examined the spectrum by fitting 
with various single component models that in- 
clude power-law (PL), blackbody(BB) and hydro- 
gen atmospheric m odel of a neutron star (NSA; 



Zavlin et al. 19961 ). We find that no single com- 
ponent can fit the spectrum well. Instead, the 
spectrum can be described by a composite model 
consisting of a non-thermal component and a ther- 
mal component. We summarize the best-fit spec- 
tral parameters in Table [TJ The non-thermal X- 



Parameter 



N H (cm"*) 

r 

F x ,pl (ergs cm -2 : 
kT (keV) 

F x ,Th (ergs cm -2 : 
R (km) 6 

X 2 /dof 



PL+BBodyrad PL+NSA C 



2.8 x 10* 1 (fixed) 2.8 x 10* 1 (fixed) 



i cn+0.10 
■ L ""— 0.09 

1.6 x icr 1;i 
0.104±0.013 
3.3 x 1CT 14 

1 QO+0.68 
1Ja -0.43 

108.9/119 



1 40+0.10 

i.6 x icr 13 
0.058t°;»;? 
3.5 x icr 14 

10 (fixed) 
108.0/119 



Note: Quoted errors indicate the 90% confidence level for 
one parameter of interest. 

a The unabsorbed flux is measured in the range of 0.5—10 
keV. 

b The radius is measured from the normalization factor 
for a source distance of 1.4 kpc. 

c The NSA model is proposed to be nonmagnetic with the 
mass and radius of the neutron star fixed at 1.4 Mq and 
10 km. 



ray emission from the pulsar PSR J0007+7303 is 
thou ght to come from t he ou ter magnetospheric 
gap (jTakata et al]|2006l [2008h . The thermal con- 
tribution with a temperature of ~ 0.6 — 1.0 x 10 6 K 
might originate from the surface of the neutron 
star. The X-ray flux from the thermal component 
is only ^20% of the total flux, which is consisten t 
with previous XMM detection (jSlane et al.ll2004l) . 
However, the observed XMM flux is ~ 4 times of 
the flux attr ibuted to the pulsar in the Chandra 
observation (jHalpern et al.l 120041) . We note that 
the aperture radius of XMM is about 15", which 
is much larger than the point-spread-function of 
Chandra and might contain emission from the 
PWN or/and the jet. The excess of the thermal 
flux detected by XMM data may be also due to 
the thermal emission contributed from the diffuse 
emission of CTA 1 dSlane et al.lll997t ). We also 



note that the Chandra spectrum was fitted with 
only 6 dofs, suggesting the model fitting might 
have large uncertainty. In this respect, the uncer- 
tainty of the thermal flux derived by the Chandra 
spectrum might be large. On the other hand, 
although the resolving power of XMM-Newton is 
worse, the larger effective area of the detectors on- 
board XMM-Newton and the longer exposures of 
the 2009 XMM data should provide much better 
photon statistics than that of the Chandra data. 
Furthermore, we note that the X-ray light curve 
shows that the pulsed components can be as large 
as ~ 20% of the DC level, indicating that the soft 
X-ray flux from the pulsar could be at least as 
large as ~ 20% of the detected flux even the ob- 
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Fig. 3. — Broad band pulsed spectrum of 
PSR J0007+7303. 
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served total flux might be contaminated by PWN 
or jets around the pulsar. We also note that the 
sum of the individual non-thermal contributions 
from the jet, PWN and the pulsar reported by 
Halpern et al.l (|2004i) is fully consistent with that 
inferred from the XMM data, though it is not pos- 
sible to disentangle various components with this 
data set. These indicate that the fluxes of these 
two data sets are not inconsistent. 

3. Discussion 

An X-ray pulsation was detected using 2009 
XMM-Newton data. However, it is difficult to ob- 
tain a pulsed spectrum with present data. The 
photons of the pulsed component are quite few 
with the pulsed fraction ~ (17±1)% of the total 
X-ray flux. In addition, the EPIC background is 
dominated by thermal emission at lower energies 
(E < 1 keV), which causes large uncertainties in 
the pulsed spectrum. For 7-ray pulsars in gen- 
eral, however, the X-ray emission can be fitted by 
a thermal component from the neutron star sur- 
face and/or a non-thermal component from the 
rclativistic particles that also emit the 7-rays. It 
has been observed from several 7-ray pulsars (e.g., 
Geminga and Vela) that the peak in X-ray light 
curve of the non-thermal component is in phase 
with the 7-ray peak, while the peak phase of the 
thermal component is different from that in the 
7-ray bands. The ~0.5 phase shift between the 
peaks of the folded light curves in the X-ray and 
the 7-ray bands, as shown in Fig. [51 indicates that 
the main pulsed component of the X-rays prob- 




Fig. 4 — Folded light curves of PSR J0007+7303 
in the 0.2-2 keV (upper panel) and in the 2-12 keV 
(lower panel). These two curves were folded with the 
frequency of 3.165844344 Hz and the frequency deriva- 
tive of —3.6133 x 1CF 12 Hz/s simultaneously at the 
same epoch zero of MJD 54897.6335129. 



ably originates from the thermal emission of the 
neutron star. 

In Fig. [31 we compare the observed pulsed 
spectrum with the theoretical predictions of 
the outer ga p accel erator model investigated by 
Takata et al.l (|2008l ). where the 7-rays and the 
X-rays are emitted via curvature radiation and 
synchrotron radiation of the relativistic particles 
in the magnetosphere, respectively. We can see 
from Fig. [3] that the observed pulsed flux in- 
tegrated from 0.2 to 12 keV of ~ (5 ± 0.2) x 



10 photons cm s is much larger than the 
predicted pulsed flux (~ 10~ 5 photons cm" 2 s _1 ) 
of the non-thermal component, suggesting the ob- 
served pulsed component is not dominated by the 
non-thermal emission from the outer gap acceler- 
ator. 

In order to examine the energy dependence of 
the pulsed X-ray emission from PSR J0007+7303, 
we divided the source photons into the soft X-ray 



■5 



(0.2 - 2 kcV) and hard X-ray (2 - 12 keV) band. 
Fig. |4] shows the folded light curves in the soft 
(upper panel) and in the hard (lower panel) X-ray 
bands. For the soft X-ray energy band, the ran- 
dom probability to yield the pulsed detection is 
3.1 x 1CT 5 [Z\ = 20.8), and the pulse profile has 
a broad sinusoidal structure, which shows a sim- 
ilar X-ray pulse profile as that in Fig. [2] In the 
hard X-ray band, on the other hand, the random 
probability to yield the pulsed detection is only 
0.039 {Z\ = 6.49), and the structure of the pulse 
profile is much more random. We did not detect 
pulsed signals at hard X-rays. This suggests that 
the X-ray pulsation detected for PSR J0007+7303 
is mainly from the thermal emission. This is also 
consistent with the fact that the thermal X-ray 
emission from a neutron star, which may not con- 
tribute significantly to the hard X-rays (> 2 keV). 
Without a firm detection of pulsed emission in the 
hard X-ray band, the nature of the power-law com- 
ponent remains unclear. 

One can use some physical properties of pul- 
sars such as ages and spin-down powers to char- 
acterize features of 7-ray pulsars. Based on these 
physical properties, PSR J0007+7303 would be 
considered to be similar to a Vela-like pulsar. On 
the other hand, one can also consider the sim- 
ilarity between two pulsars by comparing their 
high-energy emission properties. In this respect, 
PSR J0007+7303 has the certain X-ray pulsed 
detection and shares the similar spectral fea- 
tures with the Geminga pulsar; this indicates 
that the high-energy emission from these two 
pulsars may have similar mechanisms and makes 
PSR J0007+7303 as the second Geminga. Both 
the spectra of PSR J0007+7303 and Geminga 
show a steepening at the range of GcV and can 
be d escribed as a hard power- law wi t h f ~ 1.5 - 
1.6 ( Maver-Hassel wander et al. I ll994t lAbdo et al.l 
20081 ). The X-ray source spectra of these two 



pulsars can both be fitted by a composite model 
consi sting of a power-law and a blackbody comp o- 
nent (jHalpern fc Wandll997t ISlane et alj|2004) . 

However, in contrast to the Geminga pulsar, 
the X-ray spectral behavior of PSR J0007+7303 
is dominated by non-thermal contribution , whic h 



the pulsation in the soft X-ray band is expected to 
be thermal-dominant. Both sources show signifi- 
cantly dissimilar puls ed profiles for the soft X-ray 
and the 7-ray bands (jHalpern fc Rudermanlll993 
& this work). In fact, except for the Crab pulsar, 
most of the EGRET pulsars do not show resem- 
blance for the pu lse profiles in the soft X-ray and 
the 7-ray bands (jThompsonl 120011 ). The hard X- 
ray and the 7-ray pulsed profile s of the Geminga 
pulsa r have similar features (jKargaltsev et al 



might come from the PWN (plane et al] 119970 



Since the pulse profile of the pulsar centered in 
CTA 1 presents a broad si nusoidal structure a s 
that of the Geminga pulsar (jCaraveo et al. I l2004h . 



20051 ) . indicating the emission are likely originated 
from the same locations. 
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